Study of phase stability in the a-FeCr system 
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Formation energy of the cr-phase in the Fe-Cr alloy system, AE, was computed versus the occu- 
pancy changes on each of the five possible lattice sites. Its dependence on a number of Fe-atoms 
per unit cell, Nf £ , was either monotonically increasing or decreasing function of -/Vpe, depending 
on the site on which Fe-occupancy was changed. Based on the calculated AE values, the average 
formation energy, < AE >, was determined as a weighted over probabilities of different atomic 
configurations. The latter has a minimum in the concentration range where the <r-phase exists. The 
minimum in that range of composition was also found for the free energy calculated for 2000 K and 
taking only the configurational entropy into account. 

PACS numbers: 71.15.Mb, 71.15.Nc, 71.23.-k, 75.50.Bb 



Sigma (cr) phase was found in 43 binary alloys and 
many other three- or multi-component alloy systems^. 
This phase cannot be formed at the stage of solidifica- 
tion of the solution of alloying elements. Instead, it can 
only be obtained by a high temperature annealing process 
(solid state reaction). From the viewpoint of technolog- 
ical applications of alloys, the cr-phase is the one which 
should be avoided as it drastically deteriorates different 
mechanical properties of materials in which it precipi- 
tates. Sigma phase has a tetragonal unit cell (type 
P42/mnm) hosting 30 atoms distributed over five crys- 
tallographically non-equivalent sites usually called A, B, 
C, D, and E. In binary alloy systems, like Fe-Cr - the 
subject of the present paper- both alloying elements are 
present on all five lattice sites^. 

In the Fe-Cr system the cr-phase can be formed by 
an isotermal annealing in a limited range of concentra- 
tion (~45-50at% Cr) and temperature (~ 500 - 830°C)±. 
Once formed it remains stable at lower temperatures, 
but it dissolves into the a-phase at temperatures above 
830° C. Details of the reasons and mechanism of its cre- 
ation as well as that of its dissolution, hence prevention 
of its formation, are not known with sufficient clarity yet. 

Determination of the formation energy of this phase 
compared to the energy of formation of bec (a) phase 
(from which it precipitates) may allow for a better un- 
derstanding of processes involved in its formation and 
dissolution. Thus, it may be helpful in a creation of a 
new generation of materials, like stainless steels, having 
properties more suitable for a construction of new gen- 
erations of industrailly important facilities like nuclear 
power plants. 

Several theoretical papers were recently devoted to the 
a— a phase transformation^—, yet the mechanism respon- 
sible for the process is not completely understood. This 
justifies further studies toward this end. The increase of 
computing capabilities of modern computers allows cre- 
ation of increasingly sophisticated models that are able 
to include more parameters relevant to more adequately 
and precisely describe the complex structure and prop- 
erties of the cr-phase. 



In the present work we analysed semi-ordered unit cells 
having a composition in the vicinity of experimentally 
determined concentration of one of the alloying elements. 
Assumptions of the model and details of a computational 
method used in this work have been widely described 
elsewer o 9 ' 10 . 

Determination of the formation energy of the cr-phase 
has so far been based on calculations carried out for the 
unit cell which sublattices were occupied entirely by a 
single type of atom o 11 ' 12 . In the case of five sublattices 
and two kinds of atoms, there are 32 possible different 
atomic arrangements which satisfy these assumptions. 
Differences between formation energies (in fact, free en- 
talpies, H) determined for each of these configurations 
and corresponding values obtained for pure constituents 
(so called reference state - RS), AE = H a - H RS , are 
customarily presented in form of a diagram as a function 
of total concentration of alloying elements. Such pre- 
sentation enables a determination of the concentration 
range where the cr-phase exists, as well as drawing con- 
clusions regarding sublattice occupancies. Although this 
approach allows performing calculations without losing 
the symmetry of the unit cell, it does not permit to take 
into account a sublattice disorder. For example in Ref. 
12, the lowest value of the formation energy was found 
for a system named by the authors as FeCrCrFeCr. This 
notation means that sublattices A and D are fully occu- 
pied by Fe atoms, whereas B, C and E sites exclusively 
by Cr atoms. In the present paper, the convention of the 
description of the sublattice occupancy in the unit cell is 
slightly different, because we want to take into account 
atomic chemical disorder, too. Since the total number of 
atoms on each of the five sublattices is known, in our no- 
tation only the number of the Fe atoms on each sublattice 
is given explicitely, while the Cr atoms form the balance. 
Thus, according to our convention the arrangement of 
atoms FeCrCrFeCr will be called cr-20080. 

Analysis of the formation energy, as sought here, is 
fraught with a lot of problems. The most important 
seems to be the one that the cr-phase in the Fe-Cr al- 
loy system is chemically disordered, and all sublattices 
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are occupied by both alloying elements^. As found with 
neutron diffraction studies, sublattices A and D are ac- 
tually in the majority filled by Fe atoms while the other 
sublattices in the majority by Cr atoms. Yet, for the 
Fe-Cr the average sublattice occupancy is closer to the 
configuration cr-21373 rather than cr-20080. 

Another serious problem that we encounter analyzing 
the AE's is related with the fact that such analysis is per- 
formed as a function of the total concentration of one of 
the alloying elements. In fact, the AE's should be better 
considered in five dimensions i.e. versus the concentration 
for each of the five sublattices. 

Facing these problems we proposed a different ap- 
proach which so far has been successfully applied to de- 
termine different hyperfine parameters of the cr-phase in 
Fe-Cr and Fe-V systems, such as charge density, electric 
field gradients and magnetic structure of particular sub- 
lattices. It is based on the calculations performed in a 
finite number of unit cells with the symmetry reduced 
to a simple tetragonal one, in which each position (but 
not the sublattice) is occupied by one type of the atoms 
(e.g. Fe or Cr). The atoms are distributed over dif- 
ferent sites with probabilities determined from neutron 
diffraction measurements. Details of the method itself 
and the mentioned-above calculated hyperfine parame- 
ters are presented elsewher o 9 ' 10 . 

Let assume the cr-21373 be a reference state for further 
calculations. This configuration was chosen as equiva- 
lent to the average Fe-occupancy of the sublattices de- 
termined experimentally. The value of AE for the cr- 
21373 was here calculated as the average over the values 
computed for 26 unit cells of different configurations of el- 
ementary atoms on all sublattices that fulfil the assump- 
tions for the cr-21373. The next step was to calculate 
a change of AE caused by a replacement of one or two 
atoms of Fe (Cr) but on a single sublattice, only. 

The calculated in that way values of AE are shown 
in FIG la. Each point on the chart corresponds to 10- 
15 various atomic arangements analysed separately. As 
one can see, the changes of AE versus the number of Fe 
atoms, Np e , on each of the five sublattices are linear. In 
particular, increasing the amount of iron on sublattices A 
and D leads to a reduction in the AE whereas for other 
sublattices to its increase. A similar correlation could be 
observed for the cases in which the change was made on 
two sublattices at the same time without changing the 
total concentration of the alloy. An example of the latter 
is shown in FIG. lb. The results are in line with those 
presented in FIG. la, and they lead to a conclusion that 
the configuration cr-21373 is not the most energetically fa- 
vorable one for the cr-FeCr. Moreover, as it is clear from 
FIG. 1, a reduction of AE can be achieved by increasing 
the population of Fe atoms on the sublattices A and D 
or Cr atoms on sublattices B, C and E. Such changes 
of the sublattice occupancy lead to the configuration cr- 
20080, which was found as the most favorable energet- 
ically in RefJl2i. Based on the above outlined analysis 
we can conclude that if the cr-phase were created at OK 
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FIG. 1: (Online color) (a) A change of the formation en- 
ergy, AE, per one atom versus a number of Fe-atoms per unit 
cell, NFe, for different lattice sites. < AE >-values for both 
phases are shown, too. Solid lines are to guide the eye only. 
The blue band markes a concentration range of the experi- 
mentally found range of the cr-phase occurance. A different 
representation of the C- and E-sublattice data displayed for 
the er-phas in (a) is given in (b), where, additionally, an effect 
of a simultanous exchange of 2 atoms on different sublattices 
(namely, er-21274 and <r-21472) is visualized. 



(the temperature for which the calculations were carried 
out), then the sublattice occupancy would be different 
from that observed experimentally. 

One must remember that the cr-phase is a disordered 
system and each configuration set presented in FIG. 1 (as 
well as many others not taken into account in this calcu- 
lations) occurs in real terms with various probabilities. 
Therefore, in order to make a fair comparison between 
our results and the experimental ones it seems appropri- 
ate to calculate the average values of AE as a function 
of the concentration of Fe atoms. The average values, 
< AE >, were obtained based on the calculated values 
of AE for each configuration weighted by the probability 
of its occurrence. As can be clearly seen in FIG. la, the 
average AE for the cr-phase reaches its minimum almost 
in the middle of the range where the cr-phase was found 
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FIG. 2: (Online color) Magnetic, S mag , (a) and configurational, S con f, (b) entropies per one atom for the a-phase and for 
different sublattices of the a- phase versus the number of Fe-atoms per unit cell, Nfc Solid lines are only to guide the eye. 



experimentally. 

As was mentioned earlier, the cr-phase in the Fe-Cr sys- 
tem can be formed in the range of temperatures between 
~ 500 and - 830° C. Annealing above ~ 830° C leads to 
a recovery of the a-phase, whereas below ~ 500°C cr- 
phase remains stable. A comparison between the a- and 
a-phase average formation energies is presented in FIG. 
la. The < AE > a has been calculated using the KKR- 
CPA method. As one can see, its values are significantly 
lower than the values found for the cr-phase in the corre- 
sponding concentration range. This finding is consistent 
with the experimental result that the solid phase of Fe-Cr 
obtained from the liquid phase is always of the a-type. 

In order to obtain a deeper insight into a relationship 
between the energies of formation of the a- and cr-phases, 
one should not restrict itself to a comparison of the free 
enthalpy of these phases. Instead, one should also con- 
sider the free energy, F — E—TS. To this end, it becomes 
necessary to estimate the value of the entropy, S. 

The total entropy which should be taken into account 
here has several contributions, namely: configurational 
entropy, S con f, magnetic entropy, S mag , electronic en- 
tropy, S e i, and phonon entropy, S p h- 

The magnetic entropy is one of the quantities that can 
be determined separately for each of the discussed atomic 
configurations. Its contribution to the total entropy was 
calculated from the formula 

S mag =k B J2*i H^f e + 1) + (1 -Xi) ln(2/if r + 1) (1) 

where ks is Boltzmann constant, Xi stands for the con- 
centration of Fe on each sublattice, and the fii is the av- 
erage magnetic moment of Fe/Cr atoms belonging to this 
sublattice. Calculated in this work /i- values were found 
to be consistent with those already published 13,14 . Con- 
sequently, SVnag-values were obtained using eq. [T] and 
are presented in FIG. 2a. A similar for each site and 
monotonic behavior of S mag with increasing number of 
Fe-atoms, Np e , reflects the fact that the average mag- 
netic moment changes with a change of Fe-concentration 



on each sublattice in a similar way. The magnetic entropy 
calculated for the a-phase was found to be larger then the 
corresponding one of the cr-phase, but the concentration 
dependence for the two phases is quite similar. Since the 
influence of S mag on the phases formation should depend 
on the difference of the entropies, S maga — S maga , one 
should expect that this type of entropy does not influ- 
ense the concentration range of a — a transformation too 



much, because the difference S„ 



for various 



x is only weakly dependent on the concentration. 

The configurational contribution to the entropy, S con f, 
was calculated according to the known formula: 



^conf 



= -kf 



Xi hi(xi) + (1 - Xi) ln(l - Xi) (2) 



As in the case of S mag , the configurational entropy can 
also be determined here for each atomic configuration 
separately. Plots presented in FIG. 2b pertinent to the 
cr-phase reflect its complex structure. The corresponding 
values for the a-phase are larger and stay nearly constant 
in this range of concentration. Based on the results one 
should expect that these two types of entropy should have 
rather different influence on the free energy F. 

The electronic contribution to the entropy, S e i, can 
be determined from the temperature dependence of the 
electronic part of the specific heat, C v . Unfortunatelly 
the frequently used formula viz. C v = f3T + 7T 3 + ... that 
is valid as a low temperature approximation, cannot be 
used here because the temperature of ~ 500 — 800° C at 
which the cr-phase can be formed is almost twice higher 
than the Debye temperature, To, for the cr-phase (410- 
480 °C) 1 ^ and for the a-phase in the discussed range of 
concentration (~ 400°C) 1 - as determined experimentally. 
Since the high temperature changes of this type of the 
specific heat (~ /3T) are expected to be much weaker 
than other types (e.g. phonon contribution, ~ 7T 3 ) on 
one hand, and because of a saturation character of C„(T), 
on the other, we are not able to take it correctly into 
account at this stage of calculations. 
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FIG. 3: (Online color) The free energy per one atom, F, 
calculated for T = 2000 K as a function of NFe- Only the 
configurational entropy was taken into account. 



The phonon contribution, S p h, was calculated only for 
one configuration of atoms in an ordered unit cell, namely 
the one corresponding to the er-21373il. Unfortunately, 
such calculations are time-consuming and their perfor- 
mation for about 300 atomic arangements, a basis for 
this work, was in practice not feasible. The total S p h 
was also determined experimentally on the basis of a 
phonon spectrum measured on the er-phase samples for 
three different concentrations. In that case the difference 
S p h a — S p h a ~ 0.08fcs, so it is much less than that for the 
other types of entropy. Unfortunatelly, the measured val- 
ues correspond to the average entropy only, and as such 
they cannot be used for a careful analysis of the influence 
of the atomic configurations on the phase transformation. 

The above oulined analysis shows that at this stage of 
computations one has only one type of entropy, namely 
the configurational one, that can be quantitatively used 
to analyse the behavior of the free energy versus temper- 
ature. Exemplary results obtained in this study for T = 



2000 K are shown in FIG. 3. It is evident that for each 
sublattice the free energy reaches its minimum value in 
the range of concentration in which the <r-FeCr is ob- 
served experimentally. Since this phase can be formed 
only at elevated temperatures, the sublattice occupancy 
should corespond to the minimum of F at that temper- 
atures, as it is indicated by a shadowed-band. 

The temperature of 2000 K is by a factor of two larger 
than the one at which the formation of sigma occurs. 
However, the calculations were performed taking into ac- 
count the configurational entropy only, which is only one 
of several contributions to the total entropy. Taking into 
account all types of the entropy will hopefully result in 
the lowering of the temperature to the values at which 
the phase really forms. 

In summary, the results presented in this paper can be 
concluded as follows: (a) if the cr-phase could be formed 
at K the arrangement of atoms would be different (viz. 
cr-20080) than the one observed experimentally, (b) the 
average value of the formation energy at K has its min- 
imum within the composition where the cr-phase occurs, 
(c) the formation energy for the a-phase is lower than the 
one for the er-phase, and consequently always the former 
precipitates from the liquidus, and (d) among different 
possible contributions to the entropy, the configurational 
one seems to have the most significant effect on the site 
occupancy. 
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